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ABSTRACT 

We examine the ability to test the black hole no-hair theorem at the 10% level in this decade 
using the binary black hole in OJ287. In the test we constrain the value of the dimensionless 
parameter q that relates the scaled quadrupole moment and spin of the primary black hole: 
<72 = — qx 2 ■ At the present we can say that q = 1 ± 0.3 (one a), in agreement with General 
Relativity and the no-hair theorems. We demonstrate that this result can be improved if more 
observational data is found in historical plate archives for the 1959 and 1971 outbursts. We also 
show that the predicted 2015 and 2019 outbursts will be crucial in improving the accuracy of the 
test. Space-based photometry is required in 2019 July due the proximity of OJ287 to the Sun at 
the time of the outburst. The best situation would be to carry out the photometry far from the 
Earth, from quite a different vantage point, in order to avoid the influence of the nearby Sun. 
We have considered in particular the STEREO space mission which would be ideal if it has a 
continuation in 2019 or LORRI on board the New Horizons mission to Pluto. 



Subject headings: gravitation — relativity — quasars: general 
hole physics — BL Lacertae objects: individual (OJ287) 



quasars: individual (OJ287) — black 



1. Introduction 

Astronomical observations and detailed astro- 
physical considerations strongly support the exis- 
tence of black hole candidates having masses in 
the range from few M Q to few 10 10 M Q . In order 
to make sure that they actually are black holes as 
postulated in General relativity (GR), we should 
prove that at least in one case the black hole no- 
hair theorems are satisfied. 

According to the black hole no-hair theorems, 
an electrically neutral rotating black hole in GR 
is completely described by its mass M and an- 
gular momentum S (Israel 1967, 1968, Carter 
1970, Hawking 1971, 1972, see Misner, Thorne 



and Wheeler 1973 for discussions). This implies 
that the multipole moments, required to specify 
the external metric of a black hole, are fully ex- 
pressible in terms of M and S. In the case of a 
Kerr black hole, characterized by the Kerr param- 
eter x, its dimensionless quadrupole parameter q 2 
is uniquely defined by 



92 = ~X 



(1) 



where q 2 = c 4 Q 2 /G 2 M 3 and X = cS/GM 2 , 
and Q 2 is the quadrupole moment of the black 
hole (Thorne 1980; Thorne, Price and Macdonald 
1986) . 

Recently, the first attempt to probe the black 
hole no-hair theorems was made by Valtonen et al. 



1 



(2010) using the available optical observations of 
the BL Lacertae object OJ287. T he quasiperiodic 



optic al light curve of this quasar (jSillanpaa et al 



19881 ) displays temporal variations having 12 and 
60 year cycles (Figure 1). 

A simple model for OJ287 involves a secondary 
black hole orbiting a more massive primary black 
hole in an eccentric orbit having a periodic- 
ity of about 12 years, while the 60 year period 
arises from the associate d periastron precession 
( Lehto fc Valtonen 19961) . In principle, one could 



imagine many different binary models that would 
satisfy these requirements. However, there exists 
a third requirement which nails down the model. 
It is related to the observed double peak struc- 
ture in the light curve of OJ287, with the two 
peaks separated by one to two years and the pair 
occuring approximately in every 12 years. These 
observations are interpreted as being due to the 
double impact of the secondar y black hole on the 



accre tion disk of the primary (jLehto fc Valtonen 
19961) . The model was able to predict the 2007 



September 13 sharp outburst to the accuracy of 
one day (Valtonen 2007, 2008, Valtonen et al. 
2008). 

In describing the full light curve of OJ287, one 
must also calculate the indirect effects of the bi- 
nary action on the accretion flow. After the orbit 
is fully determined by the sharp impacts and the 
related short but bright outbursts, it is straight- 
forward to calculate the more gentle rise and fall of 
the light curve arising from variations in the accre- 
tion flow. Thus the rest of the optical light curve 
was also predicted with fair accuracy (Sundelius et 
al. 1997, Valtonen et al. 2009). It is a relatively 
simple matter to separate these two kinds of flux 
variations by their quite different time scales (Val- 
tonen et al. 2011, see Figure 2). 

These investigations give us the confidence to 
employ the binary black hole model of OJ287 to 
test GR. 

The first orbit model of Lehto & Valtonen 
(1996) made use of 5 outbursts, giving 4 inde- 
pendent intervals of time. They allow a unique 
solution of 4 orbital parameters, the mass of the 
primary 1.71 ± 0.15 x 10 10 M , the eccentricity of 
the orbit e = 0.678 ± 0.004, the precession rate of 
the major axis of the orbit 33.3 ± 2°, and a con- 
stant cf>o specifying the oriention of the orbit at 
some initial moment of time. 




Fig. 1. — The optical light curve of OJ287 from 
1891 to 2010. The light curve includes previously 
unpublished data obtained at Harvard by R.Hudec 
and M.Basta. The line represents the binary black 
hole model. 



pa trization flux in OJZ87 




Fig. 2. — The optical light curve of OJ287 during 
2006-2008. Only low polarization (less than 10%) 
data points are shown. There is a big "hump" last- 
ing about one year and a "spike" at 2007 Septem- 
ber 13 lasting only a few days. It is the "spikes" 
of the light curve that are used to determine the 
times of impact on the accretion disk and then to 
calculate the orbit of the secondary. 
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After the predicted 2005 outburst was observed, 
a new solution was calculated using 6 outbursts, 
allowing the determination of 5 parameters (Val- 
tonen 2007). The new additional parameter is the 
thickness of the accretion disk (scale height ~ 150 
AU), while the precession rate was updated to 
37.5° - 39.1°. The timing of the 2007 outburst to- 
gether with some new historical data allowed a so- 
lution using 9 outbursts, and solving for 8 param- 
eters (Valtonen et al. 2010). These 9 outbursts all 
follow the basic light curve shape of Figure 3, with 
a rapid rise to the maximum and then a slower de- 
cay to pre-outburst level. The time scales of the 
outbursts follow the dependence on the impact dis- 
tance established by Lehto and Valtonen (1996). 
They form a very well defined sequence. There 
are no cases when an outburst in this sequence 
was expected but was not observed. All missing 
members are at times when there were no obser- 
vations. Neither are there any extra unexplained 
members of this sequence. 

The new parameters are the spin of the primary 
black hole, with xi = 0.28 ± 0.08, the mass of 
the secondary 1.4 ± 0.1 x 10 8 M Q , and q which is 
desribed below. 

Parallel to the increase in the number of 
outbursts in the solution, the number of post- 
Newtonian (PN) terms was increased in calculat- 
ing the acceleration between the binary compo- 
nents. Valtonen et al. (2010) include the domi- 
nant order general relativistic and classical spin- 
orbit coupling, which is required in order to relate 
the dimensionless quadrupole parameter q2 of the 
primary to its Kerr parameter. They write 



'12 



-QX 



(2) 



and let q be among the 8 parameters of the so- 
lution. Its value was determined as q = 1 ± 0.3. 
Valtonen et al. (2010) noted that the timing of the 
next outburst in 2015 should help to improve the 
accuracy of the \ estimate to about ±5%. This 
conclusion along with the fact that the mass of 
the primary is determined with the accuracy of 
±1% prompted us to explore the ways of testing 
the no-hair theorem at the ±10% level in the cur- 
rent decade by measuring q more accurately. In 
the literature, there exits a number of proposals to 
test the black hole no-hair theorems, plausible in 
the next decade with the help of electromagnetic 
and gravitational wave observations. The scenar- 



2.33+4.7*(x-x0-2.59*f)/f 



Fig. 3. — The standard light curve for an OJ287 
optical outburst. It is modeled after the unpo- 
larized component of the 2007 September 13 out- 
burst. The model uses three equations in three 
sections: the formulae are written inside the fig- 
ure. The parameter / allows applications to slower 
or faster outbursts since the outburst speed is a 
function of impact distance measured from the pri- 
mary black hole. The parameter a has the value 
3. 



ios include radio timing of eccentric millisecond 
binary pulsars having an extreme Kerr black hole 
as a companion (Wex and Kopeikin 1999) and ob- 
serving several stars orbiting the massive galactic 
center black hole at milliarsec distances with in- 
frared telescopes capa ble of doin g astrometry at 
~ lO^arcseconds level (|Willll2008l ) . Further, LISA 
observations of grav itational waves from extreme 
mass ratio inspirals (jGlampedakis fc Babakl l2006) 
and quasi-normal ringdown phases associated with 
massive black hole mergers (jBerti et al.ll2006l ) will 
also try to validate black hole no-hair theorems. It 
has also been argued that the imaging of accretion 
flow around Sgr A*, if its Kerr parameter is not 
close to one, may allow the testing of the no-hair 
theor ems in the near future (jjohannsen fc Psaltisl 
201dt ). The test relies on the argument that a 
bright emission ring characterizing the flow image 
will be elliptical and asymmetric if the theorems 
are violated. Johannsen and Psaltis (2011) fur- 
ther explore the possibility of detecting modes of 
quasiperiodic variability in accretion disks as test 
cases for the no-hair theorems. 

In what follows, we briefly summarize our ap- 
proach, detailed in Valtonen et al. (2010) and list 
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the improvements desirable for the test. We then 
identify by timing experiments those impact out- 
bursts (historical and future) that will be crucial 
to constrain the q value. Finally, we discuss ob- 
servational requirements to achieve it. 

2. Additional theoretical inputs and new 
timing experiments 

8 | , 1 , 1 , , 1 
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Fig. 4. — The observations of OJ287 during the 
1983 outburst (points). The time scale is in days, 
with the zero point corresponding to 1982.964. 
The dashed line shows the standard light curve 
for an OJ287 optical outburst using / = 6. 

In order to obtain the best timing, we con- 
struct a 'standard' outburst light curve (Figure 3). 
It is based on the observations of the 2007 out- 
burst which was monitored throughout the out- 
burst both in total optical brightness and polar- 
ization (Valtonen et al. 2008). No such light curve 
was previously available. These data allow us to 
separate the underlying Bremsstrahlung emission 
(unpolarized) from the polarized flares. We model 
the underlying light curve in three separate sec- 
tions using analytical functions. From the time 
zero to time 2.59 days we use a rising power-law 
form, 

flux = 2.9(t/3) 1 - 5 r 1 - 5 , (3) 

from day 2.59 to day 3.35 we use a linear rise, 

//ux = 2.33 + 4.7(t-2.59/)/" 1 , (4) 

and beyond day 3.35 we assume a decaying power- 
law form 

flux = 5.6((t + 1.5f)/(5f)y 1 - 5 . (5) 



The flux is in mJy, and t is the time measured 
from the beginning of the outburst. 

Even though the standard light curve is adopted 
from observations, the three sections may be jus- 
tified as follows: At time zero the optical depth of 
the radiating bubble equals unity and we start see 
to the interior of the bubble. The optical depth 
decreases and larger and larger volumes of the 
bubble come to view. At day 2.59 approximately 
half of the volume is visible, and thereafter also 
the rest of the volume produces emission as fast as 
the visibility front advances into the bubble. This 
stage happens quickly, in the light travel time of 
the bubble. At the third stage the flux from the 
bubble decreases as the radiating plasma cools 
adiabatically. There is a free parameter / in the 
formulae which contracts or stretches the outburst 
time scale. This is necessary since the time scale is 
a function of the impact distance, measured from 
the primary black hole. This function is given in 
Eq. 12 and in Table 3 of Lehto & Valtonen (1996). 

To illustrate the procedure, let us take the light 
curve of the 1983 outburst. Figure 4 shows the ob- 
served points overlaid by the standard light curve. 
The value / is six, the best fitting value which is 
also in agreement with Lehto & Valtonen (1996) 
time scale. The goodness of fit is judged by mini- 
mizing the x 2 - The standard light curve is shifted 
left-to-right in order to indentify the beginning of 
the outburst, by minimizing the \ 2 °f the differ- 
ence between observations and the standard curve. 
This produces a single value for t$, the beginning 
of the outburst. In order to see how uncertain this 
value is, we have varied the observed flux values by 
±1 mJy in random uniform way. We have also var- 
ied the number of points included in the fit, start- 
ing from 40 points and adding up to 48 more points 
to the tail. In this way we get 48 values of to. 
Their distribution is centered on 1982.964 and it 
fits reasonably well with a Gaussian of a = 0.0004 

yr- 

The same procedure was applied to other eight 
outburst light curves. The results are shown in 
Table 1. Notice that the error limits are generally 
a little narrower than in Valtonen et al. (2010). 
In that paper only the first section of the standard 
light curve was used. Thus the use of the full light 
curve produces some improvement in timing. 

We may also ask how accurately the theoretical 
model can determine the outburst timing. In the 
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Fig. 5. — The disk profile during the 1912 impact 
of the secondary on the accretion disk. The disk 
central plane is shown, together with two other 
contours approximately outlining the extent of the 
disk. The data come from the particle disk simu- 
lations of Valtonen (2007). The secondary enters 
from the upper right and travels to lower left. Its 
progress is marked by ticks at the interval of 0.01 
yr. The x-axis is the radial coordinate measured 
from the primary, in astronomical units, while the 
y-coordinatc is perpendicular to the disk plane. 
The original disk was at the zero level in the y- 
coordinatc. 



impact model a typical disk crossing time of the 
secondary black hole is one week. However, events 
at much shorter time scale can make a difference 
to the outburst timing, as illustrated by Ivanov et 
al. (1998). They show that there is a factor of 
two pressure change over the distance of 1/30 of 
the disk width ahead of the secondary black hole. 
Thus we may consider time steps of ~ 1/30 of a 
week as physically meaningful, i.e. the relevant 
time step is ~ 0.0005 yr. Also in the gas bubble 
bursting out of the disk this same time scale pro- 
duces a significant amount of evolution, as shown 
in their Figure 4. We consider this the minimum 
time step that has astrophysical relevance. 

Another limitation to the timing accuracy 
comes from the influence of the secondary on the 
level of the accretion disk. The approaching sec- 
ondary lifts the disk up, and cause an impact ear- 
lier than predicted in a rigid disk model (Ivanov 
ct al. 1998). The calculation of this effect by par- 
ticle disk simulations was carried out by Valtonen 
(2007). In Figure 5 we show the profile of the 



Table 1: Outburst times with estimated uncertain- 
ties. These are starting times of the outbursts. 

1912.970 ± 0.010 

1947.282 ± 0.0005 

1957.080 ± 0.020 

1972.94 ± 0.005 

1982.964 ± 0.0005 

1984.130 ± 0.002 

1995.843 ± 0.0005 

2005.74 ± 0.005 

2007.692 ± 0.0005 



accretion disk immediately after the impact in the 
summer of 1912. We notice that the disk is lifted 
towards the approaching secondary and is bent. 
In the model we need to know the raised level of 
the accretion disk at the time of the impact. The 
orbit of the secondary is marked by ticks at the 
intervals of 0.01 yr. We estimate that the timing 
accuracy in this occasion is ±0.005 yr. This is 
the typical accuracy that we may use for impacts 
at the outer disk (i.e. the 1913, 1957, 1973, 2005, 
2015 and 2022 outbursts). For the inner disk the 
effect is negligible. For example, the 2005 outburst 
may be timed within ±0.001 yr from observations, 
but such accuracy is not justified by theoretical 
considerations. 

The PN approximation provides the equations 
of motion of a compact binary as corrections to 
the Newtonian equations of motion in powers of 
(v/c) 2 ~ GM/(c 2 R), where v, M, and R are the 
characteristic orbital velocity, the total mass, and 
the typical orbital separation of the binary, respec- 
tively. In Valtonen et al. (2010), the binary black 
hole was modeled using a spinning primary black 
hole with an accretion disk and a non-spinning 
companion. The calculation of the orbit included 
all the 2PN-accurate non-spinning finite mass con- 
tributions as well as the leading order general 
relativistic (1.5PN order) and classical spin-orbit 
(2PN order) spinning contributions, and radiation 
reaction effects (2.5PN order). Here the termi- 
nology 2PN, for example, refers to corrections to 
Newtonian dynamics in powers of (v/c) 41 . 

For the present work, we incorporated the fol- 
lowing three new features in orbit calculation. 

The non-spinning finite mass contributions to 
the binary black hole dynamics are now fully 3PN 
accurate. This is achieved by adding the 3PN con- 
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tributions to d 2 x/dt 2 , where x is the relative sepa- 
ration vector in the center-of-mass frame, as given 
by Eq. (1) in Valtonen et al. (2010b). These 3PN 
contributions are available in Mora & Will (2004). 
Secondly, we let the smaller black hole also spin. 
Therefore we add the leading order spin-spin con- 
tributions to x, given by Eqs. (54) in Barker & 
O'Connell (1979). These contributions appear at 
the 2PN order. This is a desirable addition due to 
the fact that the classical spin-orbit coupling, cru- 
cial to constraining the value of q, also enters or- 
bital dynamics at the 2PN order. Furthermore, we 
have also included the contributions due to clas- 
sical spin-orbit and general relativistic spin-spin 
interactions to the precessional equation for the 
unit spin vector si in the present analysis. The 
compact binary dynamics, employed in the present 
investigation, schematically reads 



x = 



fx 

dt 2 

dsi 
~~dt 



Xq + Xip N + x S o + x S s + x Q 

+ X 2PN + X 2 .5PN + X 3PN , (6) 

(n S o + n ss + n Q )x Sl , (7) 



The explicit expressions for the non-spinning con- 
tributions to x are listed in Mora & Will (2004). 
The spin related contributions to x and ds\/dt 
are from Barker & O'Connell (1979). The addi- 
tional spin related contributions to the dynamics, 
namely xss, ^ss and SIq, that are not listed in 
Valtonen et al. (2010a), are given by 



xss = - 



3G 3 m 3 



XiX2v\ (»i -n) s 2 



&ss 



+ (s 2 ■ n) si - 5 (si • n) (s 2 ■ n) n 
+ ( Sl -s 2 )n\, (8) 



G 2 m 2 



G 2 m 2 r] 



X2\ 3 (s 2 -n)n- s 2 > , (9) 



0) 



QXi\ 3 (»i • n)n- s{ 



where the Kerr parameter \i an( i t ne un it vector 
S\ define the spin of the primary black hole by the 
relation S\ — Gm\ xi s\/c, while Xi 1S allowed to 
take values between and 1 in GR. A similar rule 
applies to \i m ^2 = Gm 2 X2 s 2/c The vector 
n is defined to be n = x/r, where r = \x\, while 
m = mi + m 2 and r] = mi m 2 /m 2 . 



The effect of including the 3PN corrections to 
the orbital dynamics is roughly a one percent in- 
crease in the estimated mass of the bigger black 
hole, demonstrating that the employed PN dy- 
namics is in the convergent regime (Valtonen et 
al. 2010b). Therefore, its influence on the q esti- 
mate is within our desirable error limits. 



q_center=0.99, sigma=0.30 
N=1052 




Fig. 6. — The distribution of the final q- values 
among the solutions using the currently best im- 
pact timings. We combine sets 7 & 8 to create the 
g-distribution. 

The fact that the primary black hole is spinning 
slowly, in other words Xi is much less than 1, in- 
dicates that the spin-spin contributions to x enter 
the binary dynamics at the 3PN order. Note that 
it is the definition of the spin of a compact object, 
namely S <~ G m 2 v spin / c 2 , where m co and v spm 
are the typical mass and rotational velocity of the 
spinning compact object, that makes the spin-spin 
contributions to appear at the 3PN order in our 
model. Further, the presence of the symmetric 
mass ratio r/ as a common factor in the spin-spin 
corrections ( r\ ~ 10~ 3 in our binary black hole 
model) ensures that these corrections have only 
minor effects on the orbits. 

The combined effect of higher PN order and the 
presence of 77 as a common factor makes sure that 
the leading order spin-spin and classical spin-orbit 
couplings make negligible contributions to s\. The 
timing experiments also reveal that the change in 
the orientation of the secondary spin axis does not 
affect the q estimates. 

In this investigation, we make a third improve- 
ment based on astrophysical considerations. In 
Valtonen et al. (2010a), the black hole spin of the 
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primary black hole was parallel to the accretion 
disk spin at the initial epoch, which was the year 
1856. Due to PN effects, the black hole spin wan- 
ders about 9° off from this direction during its pre- 
cession cycle that lasts around 1300 years. In the 
present model, the precession cone axis coincides 
with the mean accretion disk axis. After perform- 
ing a number of numerical experiments, we found 
that it is possible to choose a suitable initial direc- 
tion for S\ such that the angle between the spin 
and the disk axes remains constant ( ~ 8°) during 
the precessional motion of si. 

It is reasonable to expect such a situation 
due to the Bardeen- Peterson effect. Because the 
time scale of the Bardeen-Peterson effect is much 
longer than the black hole spin precession time 
scale (Lodato and Pringle 2006), the two direc- 
tions do not coincide. The time scale of the 
Bardeen-Peterson effect is of the order of one mil- 
lion years (Natarajan and Pringle 1998, Eq. 2.16) 
which is intermediate between the spin preces- 
sion time scale of 10 3 yr and the binary merger 
evolution time scale of about 10 8 yr (Iwasawa et 
al. 2011). Thus we expect that in 10 8 yr the 
Bardeen-Peterson effect is important up to the 
distance of about 10 2 Schwarzschild radii in the 
disk (Natarajan and Pringle 1998, Eq, 2.8), but 
the disk can follow only the mean direction of the 
spin. It cannot keep up with the 10 3 yr evolution 
of the actual spin. 
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Fig. 7. — The correlation of xi with the zero point 
to of the 2015 outburst (labelled t(2015...). The 
insert gives the best fit functional dependence. Set 
3. 

With the above mentioned additional features, 



we have searched for orbit solutions. As before, 
an automatic search algorithm is used. It takes 
a trial orbit, then improves it until all nine out- 
bursts happen within their allotted time intervals. 
Typically one solution is found in 3 minutes of 
computing time with a modern PC. We have used 
sets of 1080 orbits with given standard parame- 
ters. However, the convergence was not always 
found in a reasonable amount of time. Then the 
attempt to find a solution was discarded and the 
next trial was started. For this reason the number 
of orbits in a set is always less than 1080. 

In Table 2 we give the set number, the orbit 
number and value of the dimensionless spin of the 
primary xi in the first three columns, respectively. 
The spin value was generally taken as xi = 0.275, 
except in two sets (3 & 11) where a range of xi val- 
ues were used. The next column in Table 2 gives 
the value of the secondary spin. The spin X2 com- 
ponents are either -0.5,-0.5,-0.5 (standard case), 
0,0,0 (set 6) or +0.5,+0.5,+0.5 (set 5). Smaller 
sets were calculated to ascertain that these three 
X2 values are representative in statistical sense of 
the different orientations and magnitudes of s 2 - 
The last column in Table 2 gives the range of 
the parameter q which is initially uniformly dis- 
tributed between the limits. The solutions con- 
verge to a distribution of q which is narrower than 
this range. Only in set 3 a fixed value of go = 1 
was used. 

Even though go is not a physical parameter but 
an ingredient of the orbit finding algorithm, its 
proper choice is still important. We tried setting 
go initially far from the value go — 1, using either 
go = or go = 2, but we found that our code 
was not able to find enough solutions to justify 
these choices. For example, in the latter case only 
23 solutions were found which concentrate around 
Qcenter = 1-16 with a standard deviation of 0.15. 
Taking the distribution uniformly between these 
two limits produces more solutions, but mostly 
from the range between go = 0.6 and go = 1.4. 
Therefore we decided to carry out most experi- 
ments using this range of go- However, since it is 
possible to add some solutions also using the wider 
range of g , we have sometimes added two sets to- 
gether, one with the narrower range, and the other 
one with the wider range. The distribution of go 
then mimics a Gaussian with the standard devia- 
tion of 0.42. The resulting q distribution is always 
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narrower than this, demonstrating that we are not 
biasing the final q distribution to be unduly com- 
pact by our choice of qo- 

There were additional conditions in some sets 
which are not listed in Table 2. In set 4 the out- 
burst uncertainty limits were taken from Valtonen 
et al. (2010). They are generally somewhat wider, 
and also some of them are centered a little differ- 
ently from the ranges listed in Table 1. (Note that 
Valtonen et al. 2010 has a misprint; one of the 
central values they use is 1995.843, not 1995.841). 
On the other hand, sets 7 and 8 explore the solu- 
tions where one of the intervals is made narrower, 
i.e. the range of the timing of the 2005 outburst 
is 2005.74 ± 0.001, five times norrower than in our 
standard case. In general the initial angle between 
the disk and the primary spin xi is 8°, but in set 
10 we tested also the case of an initially zero angle. 
In sets 11 and 12 the to in 1995 is shifted down and 
up by 3.5 hours, respectively. 

3. Results 
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q=3.6-718*(t-1959.21| +/-0.018 



Fig. 8. — The correlation between q and the 
starting time of the 1959 outburst t (labelled 
i(1959.21..) whcnxi = 0.275. The functional form 
of the correlation is given inside the figure. There 
is also a correlation with the \\ such that the line 
of regression is shifted to the right in the figure by 
one unit for a decrease of xi by 0.015 units. Five 
units in the time axis corresponds to 4.4 hours. 



For every set we have constructed the distri- 
bution of q values, and since these distributions 
resemble a Gaussian, we have determined the 
best fitting Gaussian parameters, the central value 
Qcenter and the standard deviation a for each dis- 
tribution. These are listed in Table 3, together 
with the errors in each parameter. Figures 6 illus- 
trates one such distribution, a combination of sets 
7&8. 



Table 2: Parameters of the 



Set 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 



No. 

1012 

864 

901 

362 

1009 

1017 

598 

454 

914 

283 

658 

1015 



Xi 
0.275 
0.275 

0.26±0.04 
0.275 
0.275 
0.275 
0.275 
0.275 

0.26±0.03 
0.27 
0.275 
0.275 



sets of 

X2 

-0.87 
-0.87 
-0.87 
-0.87 
+0.87 
0.0 
-0.87 
-0.87 
-0.87 
-0.87 
-0.87 
-0.87 



solutions. 
?o 

0.6-1.4 
0.0-2.0 
1.0 

0.6-1.4 
0.6-1.4 
0.6-1.4 
0.6-1.4 
0.0-2.0 
0.6-1.4 
0.6-1.4 
0.6-1.4 
0.6-1.4 



The Kerr parameter of the primary black hole 
Xi will be constrained by the timing of the out- 



burst in 2015 (Valtonen et al. 2010a). Figure 7 
shows the correlation of \i with to, the zero point 
of the 2015 outburst, using set 3. The accuracy of 
the xi determination, after the the 2015 outburst 
time is known, is ±0.005 (one cr). It is likely that 
OJ287 becomes a "Christmas star" of 2015. 

We will now discuss those historical outbursts 
which were not employed in finding the orbital so- 
lution. We will look for a correlation between the 
starting times of these outbursts with q. We will 
also ask whether the distribution of q can be nar- 
rowed down by future observations. 

The first outburst which we studied is the 
1934 outburst, with the expected starting time 
1934.3439 if X i = 0.275 and q = 1. The corre- 
lation of the starting time with q is so weak that 
this outburst is of no interest in determining q. 
Moreover, there are no data in the historical light 
curve yet to verify this outburst. The 1935 out- 
burst is no better in this respect. It is expected 
at 1935.3939, with little correlation with q. This 
outburst has not been verified either in the obser- 
vational record. 

The next interesting outburst should have 
taken place in 1959. Here we expect quite a strong 
correlation between the start of the outburst and 
q (Figure 8). 
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Tabic 3: Gaussian fits to q distributions. 



bet 


Qcenter 


error 


(J 


error 


1 
1 


l.UU 


U.Ul 


n op: 
U.ZO 


U.Ul 


z 


1 m 
1.U1 


U.UO 


u.oy 


U.U4 


Q 

6 


1 HQ 


U.Ul 


U.Uo 


U.Ul 
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5 


0.96 


0.01 


0.24 


0.01 


G 


0.98 


0.01 


0.25 


0.01 


7 


0.98 


0.01 


0.25 


0.01 


8 


0.98 


0.02 


0.42 


0.02 


9 


1.00 


0.01 


0.27 


0.01 


10 


0.58 


0.01 


0.19 


0.01 


11 


0.74 


0.01 


0.19 


0.01 


12 


1.06 


0.01 


0.30 


0.01 


1+2 


1.02 


0.01 


0.33 


0.01 


1+2+5-9 


0.99 


0.01 


0.12 


0.01 



We see from Figure 8 that the difference in the 
starting time of the outburst by 15 units corre- 
sponds to the range of 0.8 in q. Since 15 units in 
the figure corresponds to 13 hours, with good tim- 
ing, say with 5 hour accuracy, it should be possible 
to determine the q value at the level of 1 ± 0.20 if a 
fair number of detections are found in the histor- 
ical plate collections. The current observational 
situation is depicted in Figure 9. 

The next outburst of interest is the 1971 out- 
burst and at present there exists only one observ- 
ing point (Figure 10). Figure 11 shows the ex- 
pected correlation between q and to- If enough 
observations are found, it should be possible to 
determine the q value with the accuracy as high 
as 1 ± 0.16 units. 

The 1995 outburst was already used in our so- 
lution. There was an intensive monitoring cam- 
paign of OJ287 (called OJ94) during this outburst 
season, but unfortunately there exists a gap in 
these observations just at the crucial time (Figure 
12). It may still be possible that there are mea- 
surements somewhere which are not recorded in 
the OJ94 campaign light curve, and which would 
be valuable in narrowing down q even from these 
data. The line in Figure 12 is drawn using the 
standard light curve of Figure 3 as a template to 
compare with the 1995 observations. In set 11 the 
value of tg has been shifted down by 3.5 hours, and 
in set 12 it is shifted up by the same amount. The 
shifts lead to shifting q cen ter up or down by ~ 15%. 
It should be noted that even a few more measure- 




Fig. 9. — The observation of the brightness of 
OJ287 at the expected 1959 outburst time. The 
dashed line is the template from the 2007 outburst 
while the squares reprent the three observations. 
The base level has been arbitrarily normalised as 
there are not enough data to determine it. 

mcnts of 1995 could narrow down the range of q. 

Let us now turn our attention to the expected 
future outbursts in our binary black hole model. 
We expect three more outbursts during the next 
two decades occuring in 2015, 2019 and 2022. As 
we mentioned above, the 2015 outburst should be 
an easy one to detect, as it is expected in Decem- 
ber of that year. The exact date will in fact give 
us a good spin value. The dependence on q is sec- 
ondary, and thus it is of no use by itself for the 
testing of the no-hair theorems. 

The 2019 outburst is sensitive to the q and with 
X = 0.275, it should begin at 2019.53175 if q = 1 
(Figure 13). With good timing the q value is de- 
termined with the accuracy of 1 + 0.16 (Figure 14). 
If by good luck we will find the necessary histori- 
cal data to time both the 1959 and 1971 outbursts 
in addition to observing the 2019 outburst, we get 
close to the 10% accuracy in q (see the combined 
sets 1+2+5-9 in Table 3). Improvements in under- 
standing the astrophysical processes in O J287 may 
then bring the accuracy even below 10%. Without 
the accurate timing in 2019, the q value cannot be 
determined better than to 1 ± 0.3 even if the Xi 
determination in 2015 is a success (Figure 6). 

However, observing OJ287 during the expected 
2019 outburst window, namely around 2019 July 
21, is practically impossible from the ground. This 
is because the angular distance between the Sun 
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Fig. 10. — Observations of the brightness of O J287 
at the expected 1971 outburst time. The solid 
line is the template from the 2007 outburst while 
the square represents the sole observational point. 
The base level (dashed line) is uncertain. 



and OJ287 in the sky is only about 12 degrees at 
the beginning of this event, and it goes down to 8 
degrees by the time of the peak flux. It may re- 
quire space observations to carry out the measure- 
ment. The 2022 outburst is scheduled practically 
at the same time of year as the 2019 outburst. 
Obviously it would also be of interest to observe 
this event as it would tie down the parameters of 
the general model more narrowly. However, it will 
not give any further information on q. The 1922 
outburst follows an impact on the outer disk, and 
these impact timings are not sensitive to q. 

Finally, let us comment on the effect of addi- 
tional features introduced in the present study. 
First, the improvement in the timing of the out- 
bursts with respect to Valtoncn ct al. (2010a) 
does improve the accuracy of the q determina- 
tion. Comparing set 4 with sets 1, 5 and 6 we see 
that the a of the distribution is greater by ~ 10% 
when the timing intervals are wider. If the angle 
between the disk axis and the primary spin axis 
varies as in Valtonen (2010a), the distribution is 
not centered on q = 1 (set 10), unlike in our stan- 
dard case. The solutions in Valtonen (2010a) were 
not numerous enough to detect this effect reliably. 

The effect of the \2 on the q distribution was 
calculated in three cases: first with the "normal" 
direction of the secondary spin \2 , then with the 
opposite spin, and finally with zero spin. The 
three distributions are different, showing that the 
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q=4.9-503'(t-1971.12)+/-0.14 



70 75 80 

1(1971.12..) 



Fig. 11. — The correlation between q and the 
starting time of the 1971 outburst t (labelled 
t(1971.12..) whcnxi = 0.275. The functional form 
of the correlation is given inside the figure. There 
is also a correlation with the xi such that the line 
of regression is shifted to the right in the figure by 
one unit for a decrease of xi by 0.01 units. Five 
units in the time axis corresponds to 4.4 hours. 



spin-spin interaction has influenced the orbits (sets 
1, 5 and 6), but the Gaussian parameters of the 
distributions are only marginally different from 
each other. 

4. The methods of observation of stars 
near to Sun and the possibilities to ob- 
serve OJ287 in 2019 

The objects at small angular distances from 
the Sun (the estimated value for OJ287 is 8 to 
12 degrees) are difficult to observe due to high 
background caused by intense sunlight. The vi- 
gnetting of the highly luminous solar disk enables 
a reduction of the background and observion of 
stars near the Sun. Recently two different meth- 
ods have been used to observe stars at small an- 
gular distances from the Sun, namely the corona- 
graph method, and the hclioscopic imager method. 

4.1. Coronagraph method 

The coronagraph has been in use already for a 
long time both in the ground based as well as in 
space based observations. The SOHO's LASCO 
C3 can serve as an example of a recent corona- 
graph in space (Morrill et al. 2006). The SOHO 
spacecraft has three coronagraphs (LASCO) on- 
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Beginning of the 1995 outburst 
at 1995.843 



Day -Julian Day 2450026.65 



Fig. 12.— The observations of OJ287 at the be- 
ginning of November 1995, transformed to opti- 
cal V-band. Altogether 50 observations have been 
binned to 7 points. Overlaid is the theoretical light 
curve profile from Figure 3. The zero point of time 
is at Julian Day (JD) 2450026.65, i.e at 3:36 hours 
GMT on November 6, 1995. 



board, two of which are still working (C2 and 
C3) . We have used the publicly available images of 
these experiments to estimate the expected limit- 
ing magnitudes for this method. Estimating lim- 
iting magnitude for C2 is difficult because most 
of the LASCO C2 image is obscured by the solar 
corona and only a few stars are visible. LASCO 
C3 covers an area of 32 diameters of the Sun (i.e. 
about 16 degrees, hence the OJ287 position dur- 
ing the 2019 predicted flare would be covered) and 
stars are clearly visible in the images, although a 
large part of the image is also obscured by the 
solar corona. 

There are several problems which complicate 
the estimation of the limiting magnitude. The 
dominant one is the stray light which can be mis- 
taken as background stars. There also several 
problems in determining the position and rotation 
of the spacecraft. One has to align the stars on 
the image from LASCO C3 and a star chart for 
the position (AAVSO charts can be used). For 
the elimination of the stay light, the video which 
is provided on the SOHO web pages was used. If 
there is a point source in five images in a row, we 
consider it a star. The rotation also complicates 
determination of the background stars. The best 
method here is to find a noticeably bright star and 
to align the image with the chart containing this 
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q=1. 76-438*(t-2019.53) +/-0.14 
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Fig. 13. — The correlation between q and the 
starting time of the 2019 outburst t (labelled 
£(2019.53..) when xi = 0.275. The functional form 
of the correlation is given inside the figure. There 
is also a correlation with the xi such that the line 
of regression is shifted to the right in the figure by 
one unit for a decrease of xi by 0.0125 units. Five 
units in the time axis corresponds to 4.4 hours. 



star and its surroundings. 

Pleiades arc possibly the best objects to use, be- 
cause the star cluster clearly defines the rotation. 
For our work we have used an image of Pleiades 
where stars show up to the limiting magnitude of 
10. It is compared with an the image from LASCO 
C3 (14.05.2010, Sun approaching Pleiades). The 
faintest stars of Pleiades detectable in the LASCO 
C3 image are of magnitude 8. Hence at best, we 
can obtain a limiting magnitude around 8 at the 
edges of the field where the corona is faint. Also 
the result strongly depends on the state of corona 
as during strong coronal mass ejections the limit- 
ing magnitude will be lower. 

Probably a deeper magnitude could be achieved 
for a space based coronagraph with a larger aper- 
ture, but this would require an independent fea- 
sibility study as the previous coronagraphs were 
designed for solar studies, not for photometry of 
nearby stars. 

In an independent study of the C3 limiting mag- 
nitude by Andrews (2000) using the same target 
(Pleiades cluster) , deeper limits of magnitudes be- 
tween 10 and 14 were achieved. As we would re- 
quire a limiting magnitude of 15 or better for our 
timing measurement in July 2019, it is clear that 
a LASCO C3 type instrument will not be able to 
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0.4 0.6 0.8 1 1.2 1.4 1.6 

q 

Fig. 14. — The distribution of q in the simulations 
where the 2019 outbursts has been timed with the 
accuracy of 5 hr and also the %i value has been 
found by the observations of the 2015 outburst. 
This gives the practical limit of the method of us- 
ing OJ287 to test General Relativity if historical 
records do not bring additional information at key 
moments of time. 

do the job. 

4.2. The imager method 

Another method recently used in a space exper- 
iment is the method of heliospheric imager which 
is part of the SECCHI experiment (Howard et al. 
2008, Eyles et al. 2009). The SECCHI experi- 
ment is onboard the STEREO space mission and 
consists of five telescopes, which together image 
the solar corona from the solar disk to beyond 1 
AU. These telescopes are: an extreme ultraviolet 
imager (EUVI: 1-1.7 solar radii), two traditional 
Lyot coronagraphs (COR1: 1.5-4 solar radii and 
COR2: 2.5 - 15 solar radii) and two new designs of 
heliospheric imagers (HI-1: 15-84 solar radii and 
HI-2: 66 - 318 solar radii). All the instruments 
use 2048 x 2048 pixel CCD arrays in a backside- 
in mode. The EUVI backside surface has been 
specially processed for EUV sensitivity, while the 
others have an anti-reflection coating applied. 

The HI objectives, like the rest of the SECCHI 
suite, make visible light observations of CMEs and 
other structures as they transit from the corona 
and into the heliosphere. The HI package consists 
of two small, wide-angle telescope systems (HI-1 
and HI-2) mounted on the side of each STEREO 
spacecraft, which together view the region be- 



tween the Sun and the Earth. HI has no shut- 
ter mechanism, other than a one-shot door that 
protects the instrument from contamination dur- 
ing ground operations and the launch. Thus, an 
image is collected in a shutterless mode, in which 
the intensity at each pixel is an accumulation of 
the static scene and a smearing of the image dur- 
ing readout. This smearing can be removed on the 
ground. 

The HI instrument concept was derived from 
the laboratory measurements of Buffington et al. 
(1996) who determined the scattering rejection as 
a function of the number of occulters and the an- 
gle below the occulting edge. The concept is not 
unlike observing the night sky after the Sun has 
gone below the horizon. 

While the specific concept used here has not 
been flown before, two other instruments have 
flown which have validated the ability to mea- 
sure the electron scattered component against the 
strong zodiacal light and stellar background. The 
Zodiacal Light Photometer (Pitz et al. 1976) on 
the Helios spacecraft, launched in 1974, and the 
Solar Mass Ejection Imager (SMEI) instrument 
(Eyles et al. 2003), on the Coriolis spacecraft, 
launched in 2003 have demonstrated that a prop- 
erly baffled instrument can detect CMEs (Tappin 
et al. 2003). 

The HI-1 and HI-2 telescopes are directed to 
angles of about 13 degrees and 53 degrees from 
the principal axis of the instrument, which in turn 
is tilted upwards by 0.33 degrees to ensure that 
the Sun is sufficiently below the baffle horizon. 

The novel heliospheric imagers achieve magni- 
tude limits for stars of about 13 to 14 in a 40 
minute exposure. The HI-1 imager covers the re- 
gion of 7,5 to 24 degrees from Sun and is hence well 
suited to observe the OJ287 during the predicted 
2019 event (OJ287 is expected to brighten from 
~ 14.3 to ~ 13 in V-magnitudc). A specific design 
for a dedicated space experiment to observe OJ287 
optimizing the performace for stellar fotometry at 
angular distances around 10 degrees from the Sun 
may lead to a limiting magnitude increase to about 
15. 

Another possibility involves pointing Long 
Range Reconnaissance Imager (LORRI) in the 
'New Horizons' mission to Pluto that consists of 
a telescope with a 20.8 centimeter aperture at 
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OJ287 for one week in July 2019. In the 4 x 4 
pixel binning mode the limiting magnitude in V 
is expected to be greater than 17 (Cheng et al. 
2008). 

4.3. Secondary science 

The OJ287 heliosphcric imager is expected to 
provide in addition to the OJ287 photometry valu- 
able data for other scientific fields such as moni- 
toring of astrophysical targets near to Sun, and 
optical searches for optical counterparts to GRBs 
occurring at small angular distances from the Sun. 

5. Discussions and conclusions 

We have shown that it is possible to test GR 
at the second PN order using the binary black 
hole system in OJ287. We find that GR can be 
confirmed with the one sigma accuracy of 30% us- 
ing the currect observations and theoretical un- 
derstanding of the system. One of the theoretical 
conditions is that the rotation axis of the accre- 
tion disk is at a constant angle with respect to the 
precessing spin axis of the primary black hole. It 
may be possible to verify this in future by study- 
ing the structure of the radio/X-ray jet in OJ287, 
and by theoretical studies of how the jet direction 
is determined when the two axes are not parallel 
to each other. With this proviso, we argue that it 
should be possible to test, in principle, the black 
hole no-hair theorem at the 10% level in the cur- 
rent decade by employing the binary black hole 
model of OJ287. 

We have also shown that the third order PN 
terms are too small to be detected in the OJ287 
system. They depend on the exact value of the pri- 
mary mass, at the 1% level, as well as on the spin 
of the secondary, for which it is difficult to find in- 
dependent measurements at the required level of 
accuracy. 

We demonstrate that the testing at the above 
precision will require a certain amount of good 
luck in the sense that there should exist some yet 
unknown observations in the historical records at 
certain crucial time windows. Also it is highly 
desirable to have space-based optical observations 
to monitor the impact outburst of OJ287 in 2019. 
One possibility is to use the 'New Horizons' mis- 
sion to Pluto which by 2019 is already past Pluto. 

Employing a suitable solar observing mission 



to monitor OJ287 could be another option. We 
have looked at the SOHO coronagraph images and 
find that the limiting magnitude there is ~ 8. We 
would need to get to magnitudes <~ 15 in blue or 
UV, and thus one would need to cover the inner- 
most 5 degrees of the solar image instead of just 
the solar disk. A better case for the no-hair test 
would be a continuation or follow-up mission of 
STEREO which would for one week concentrate 
on OJ287 instead of solar flares. 

Due to the fundamental nature of the test and 
the fact that the astrophysical systems associated 
with the other proposed tests are yet to be ob- 
served and likely to be plausible only in the next 
decade, it may not be even be out of question to 
plan a small space mission to monitor the 2019 
outburst and hence to test the black hole no-hair 
theorems. 

R.H. acknowledges the grants 102/09/0997 and 
205/08/1207 by the GA CR and ME09027 by 
MSMT. 
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